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■ Relativistic integral representation in terms of experimental neutron- 

■ proton scattering phase shifts is used to compute the charge form fac- 



tor of the deuteron GcdiQ 2 )- The results of numerical calculations of 
\Gcd(Q 2 )\ are presented in the interval of the four-momentum transfers 
squared < Q 2 < 35 fm~ 2 . Zero and the prominent secondary maximum in 
O ■ \Gcd(Q 2 )\ ar e the direct consequences of the change of sign in the experi- 

mental 3 Si- phase shifts. Till the point Q 2 ~ 20 frrT 2 the calculated total 
relativistic correction to \Gcd(Q 2 )\ is positive and reaches the maximal value 
of 25% at Q 2 ~ 14/m" 2 . 
Deuteron is the brightest example of intersection of nuclear and particle physics. During 
more than sixty years it serves as a source of important information about the nuclear forces, 
mesonic and baryonic degrees of freedoms in nuclei, relativistic effects and a possible role of 
quarks in nuclear structure. Therefore it is not surprising that currently the electromagnetic 
(EM) structure of the deuteron is a subject of intensive theoretical (the list of publication 
is immense) and experimental investigations. 
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With new experimental data from Jefferson Lab on elastic electron-deuteron scattering 
I ]J, |0 at momentum transfers in the GeV-range, one needs to develop relativistic approaches 
to the (np)-bound state problem. Recent experimental results from MIT-Bates provided 
the first experimental evidence for a zero in the deuteron charge form factor Gcd a t about 
Q 2 = 20 fm -2 predicted in a number of theoretical models (or not predicted, as in some kinds 
of quark models). 

Here we report the new results of numerical calculations of relativistic deuteron charge 
form factor (RDCFF). These calculations are based on the approach to the relativistic 
impulse approximation, which was briefly discussed in ref. Q (see also the review |j] and, 
especially, the references therein). The more detailed formulae are contained in ref. ||. In 
this approach the deuteron electromagnetic form factors (EMFF) are expressed in terms of 
experimental neutron-proton (n — p) phase shifts in the triplet channel of n — p scattering 
and nucleon EMFFs. The approach was developed during a number of years and is still 
being developed. The key components are as follows. 

1. General relativistic parametrization of the matrix elements of the EM current between 
the n-p scattering states (n'p'\j^\np) in terms of the invariant EMFFs of the n-p system. 
This parametrization is carried out in the canonical basis \P [iJlSmj) of the two-particle 
state vectors. Here = (p^ + n M ) is the total 4-momentum of n-p system; J, mj are the 
total angular moment and its projection; I, S are the orbital moment and total spin of n-p 
system. As the result for the n-p scattering state in the elastic channel with deuteron's 
quantum numbers (J = S = J' = S' = = 0, 2) we obtain sixteen invariant FFs 

Gf' (s, s' ,t, ...), where s = {p^ + n^) 2 , t = q 2 = — Q 2 < 0Q and the index i numerates the 
charge(C), magnetic(M), quadrupole(Q) and the quadrupole of the second generation, or 
toroidal(T) FFs. 



For the choice of kinematical variables see App.A. 
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2. The decomposition of the total G into the sum 

G = g + G int , (1) 

where g" is the FFs of unconnected part of the current (see Fig.l). This part is defined as 
usual: 

WiM\p)(n'\n) + {n'\j^\n){p'\p). 

The second term in eq.(|l|) describes the "real" n-p interaction. The decomposition in eq.(|IJ) 
is relativistic invariant due to the existence of the one-particle singular invariants of the type 
2p 8(p-p'). 

3. The calculation of g ~ as a function of invariant variables and nucleon EMFFs. These cal- 
culations were done by the methods of relativistic kinematics (without any v /c-expansion) . 
It's also evident from the two previous points that the nucleon EMFFs,entering to g ", de- 
scribe the nucleons on their mass shell. 

4. The local analyticity of the FFs 

G mt = Bm G* nt (s ± ie, s' T iv, •••) (2) 

in the vicinity of the physical region of the variables s, s' > AM 2 . As we defined above, 
s and s' are the squared invariant masses of n-p system in initial and final states. 

5. The representation of the n-p scattering matrix S(s) in the physical region in terms of 
the Jost matrix B(s): 

S(s) = B(s-ie)B- 1 (s + ie), s > AM 2 , (3) 

i.e. the solution of the boundary problem 

S(s)B + (s) = B.(s), s>4M 2 (4) 

in the scattering theory. As it well known, in the model one-channel case the the solution 
of eqs.(|3],[4D has the simple form 
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K J \ S- AM 2 J 1 \ 71 J4AP s-u J 

where 5 is the phase-shift for one channel ^-matrix (S = exp(2i<5)). But in the channel of 
n-p scattering with the deuteron's quantum numbers the scattering matrix is 2 x 2 matrix: 

t cos2e • e 2iS isin2e • e i(s+r,) ^ 



S(s) = S[5,rj,e] 



\isin2e • e i(<5+r?) cos2e-e 2i ' 7 j 



(5) 



In this case the solution of the boundary problem (eq.(^j)) is much more complicated, but it 
appears to be possible to find B in the form of series with fast convergence (see below). It 
is essential that every element of the inverse matrix B~ x has a simple pole at the point of 
bound state s = Mj. 

6. Solution of the linear conjugation problem for G \ Relative to variable s, for example, 
the initial equation is 

[g(s, ..) + G(s + it, ..)] = [g(s, ..) + G(s - ie, ..)] • S(s). (6) 

Following by this way and combining the solution of eqs.(|3|, ||) for s- variable and the analo- 
gous one for s'- variable, we can obtain the integral representation for the relativistic EMFFs 
of n-p system in terms of g and B. The resulting formulae are cumbersome enough and we 
omit them here. 

7. Taking the residues at the deuteron poles s = s' = Mj we obtain the final formulae for 
the DEMFFs (see,for example, eqs.(|7|,|8|). 

After this sketch of the method we return to the analysis of RDCFF. Among the three 
DEMFFs the DCFF is the most interesting for us when investigating the deuteron struc- 
ture. The point is that the most informative features of the deuteron structure — zero (or 
" diffractional minimum" ) and prominent secondary maximum in the theoretical predictions 
for DEMFFs — appear already in the DCFF at not so large values of Q 2 (~ 20 — 30 fm~ 2 ). 
The measurements of \G c d(Q 2 )\ for such Q 2 are real [1,2]. On the contrary, the theoreti- 
cal picture of DM(Q)FFs is so that indicated above the fine structure in this FFs appears 



for relatively high Q 2 , which hardly can be reached in current and planned polarization 
experiments. 

Following || the formulae for RDCFF G c d{Q 2 ) appear as 

-( P 5 20 + 5 22 )(p5°° + B^)[ G ^m + c™m\ + 

HPB 00 + B 02 fG? d (Q 2 ) , (7) 



G 



oo 



dsAB\s) M S D ds'g c (s,s't)AB(s') 



a' _ p2 

cd * am* s-Mj J Sl (s,t) s'-Mj 



s 2>1 = 2M 2 + ^(2M 2 - t) ■ (s - 2M 2 ) ± 



-t)(4M 2 -t)s(s-AM 2 ) . (8) 



2M 2 

In eq.( 0) p is the constant, which describes the mixing of two n — p states with different 
orbital moments (1 = and I = 2) in the point of the bound state, i.e., the deuteron. This 
constant is defined by the correspondence principle. Analyzing the nonrelativistic limits of 
eqs.(^,||), we can prove that p appears to be the standard asymptotic D/S -ratio of the radial 
deuteron wave functions, so p = 0.0277 (numerical calculations show that the dependence of 
DCFF on the variation of p is very weak). All four elements of the matrix B ll '(s) (I, I' = 0, 2) 
are taken at the bound state point s = Mj (M<j = 2M — e, where Ma, M being deuteron 
and nucleon masses, and e- being the deuteron binding energy). 

In eq.(^|) T 2 is the normalization constant, which is calculated from the condition 
G cd {0) = 1. Matrix functions AB ll '(s) = B u '(s + ie) - B 11 ' (s — ie) are the discontinuities of 
the Jost matrix B(s). The reduced Jost matrix B in eq.( [?[) is the solution of the same as for 
B, eq.( H) with the scattering matrix S = S 5,e,fj . The reduced phase shifts 5,e, f] have 
the auxiliary nature and describe the n-p scattering without bound state (S(E — >• 0) — >• 0). 
Roughly speaking, matrix B is the rest of matrix B after taking the residue at the point of 
the deuteron's pole. The formulae for B and B in terms of n—p phase shifts are cumbersome 



and are summarized in Appendix B. All relativistic aspects of the two-nucleon problem are 
contained in G w - matrix. 

As was mentioned above, the matrix functions g l f (s , s' , t) of three variables are the 
relativistic CFF of the unconnected part of the matrix element of EM current {n'p'lj^np). 
The results of the calculations of g l J:' are given in Appendix C. It is interesting to note that 
in general case in the relativistic region g l * — functions are not factorized in s, s' variables, 
whereas in the nonrelativistic limit such factorization takes place. It means that in the 
framework of the used relativistic approach || it is impossible to introduce some kind 
of concept of relativistic deuteron wave function. 

The n — p phase shifts were taken from the recent analysis of Virginia Tech group ]7[ 
and is shown in Fig. Q This analysis was made in the energy range < Ei a t < 1300 MeV. 
Extrapolation to higher energies is not so important for the calculations of Gcd for the 
small and intermediate values of Q 2 (see below). The only essential circumstance is that 
3 5i-phase shifts change sign from positive to negative and have the minimum near the 
energy E^ ~ lGeV, then go to zero in accordance with the Levinson's theorem. Any 
realistic n — p 3 Si - phase shift analysis has such a behavior. Two other states ( 3 Di and 3 Ei 
) give a relatively small contribution to Gcd- 

For the calculations of Gcd we used (as a first step) the simplest choice of the nucleon 
form factors: Ge p = (1 + Q 2 /18.23 fm~ 2 )~ 2 ,GM P / = Gun/^n = Ge p , Geu = for all 

Q 2 - 

The result of the calculations are presented in Fig. [| Our brief conclusions are as follows. 
1. The appearence of zero and secondary maximum in \Gcd(Q 2 )\ at intermediate values 
of Q 2 is the direct consequence of the change of sign of the experimental 3 S*i- phase shifts 
at intermediate energies. It is easy to calculate that the model's S(E), which decreases 
monotonically with E and is always positive {6(E) > for all E), immediately leads to 
monotonically decreasing with Q 2 values of \Gcd{Q 2 )\ without any fine structure. 
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2. Almost up to the point of zero (Q 2 ~ 20 fm~ 2 ) of \Gcd(Q 2 )\ the total relativistic 
correction (TRC), i.e., the difference between Gcd calculated relativistically (1,2) and its 
nonrelativistic limit, is positive and appears to be not small. For example, for Q 2 ~ 14 fm~ 2 
it reaches the value of 25%. 

3. TRC becomes large in the region of the secondary maximum of \Gcd{Q 2 )\, increasing the 
magnitude of the form factor. 

4. The obtained results are consistent with the available data on Gcd from MIT-Bates [[|. 
New data from Jefferson Lab E-94-018 |T[ are extremely important to test the proposed rel- 
ativistic approach in the region of higher transferred momenta, where relativistic corrections 
appear to be significant. 

We would like to make the following comments to the obtained results. Firstly, as it 
was mentioned above, for the fixed till Ei ab = 1300 MeV 3 Si -phase shifts set of ref. J7| the 
character of extrapolation to high energies is not essential when calculating the RDCFF. 
Two other, than in Fig.^| possible extrapolations of 5(E — > oo) are shown in Fig.f|. Their 
foundations are the following. If we include the so-called forbidden bound states in NN 
system (possible in some kinds of theories) in the formulation of Levinson's theorem, then 
the immediate consequence for 5(E) would be 5(E — > oo) — ► — 7r (see, for example, ref. 
||. The second extrapolation is connected with the phase shift analysis of proton-proton 
scattering for the several Ge V region (ref. [0] ) . It is natural to think that for high energies the 
contribution of Coulomb effects are small, so n-p and p-p phase shifts in the equal spin states 
are very close to each other. The results of the calculations of \Gcd(Q 2 ) I with the phase shifts 
from Fig.^| are shown in Fig.^. It is evident that the difference between these three curves 
is negligible. But when we turn to other phase shifts, which are differing from the results 
of ref. @ in the intermediate energy region, E^ < 1 GeV, the situation change essentially. 
The dependence of \Gcd(Q 2 )\ structure on the choice of different sets of experimental n—p 
phase shifts available from the literature (ref. ||10|| ) is strong enough. Possible variation of 
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5, e,rj may shift the position of zero in \Gcd(Q 2 )\ from the indicated point Q 2 = 20 fm~ 2 
to the point Q 2 = 16 fm~ 2 or to the point Q 2 = 23 fvrC 2 . At the same time the secondary 
maximum is located in the interval 26 < Q 2 < 32 fm~ 2 , and its height may change by a 
factor of seven. We can see that for improving our understanding of \Gcd(Q 2 )\ it would be 
desirable to have a more detailed phase shifts analysis of n — p scattering in triplet channel 
in intermediate energy region E iab < 1 GeV. Secondly, let us indicate the dependence of 
\Gcd{Q 2 )\ on the possible choice of nucleon EM form factors. Since the uncertainties of 
Ge p (Q 2 ) in the considered range of Q 2 are very small, the main effect in \Gcd(Q 2 )\ m &y be 
caused only by variation of GEn{Q 2 )- It seems to be generally accepted that the maximal 
deviation of Geti{Q 2 ) from the zero- value approximation Geu = is given by known formula 
Geu{Q 2 ) — — V-n T G Ep{Q 2 ) j where (i n = —1.91 is the neutron anomalous magnetic moment 
and r = Q 2 /AM 2 . The results of the calculations of \Gcd{Q 2 )\ with this nonzero values of 
Geu(Q 2 ) are shown in fig. 2. One can see that the effect is sizable and the contributions of 
relativistic effects and nonzero Geu have a similar behavior. 

Finally, we show for comparison in Fig.|3| the results of calculation of Gcd i n a relativistic 
approach developed in ref. 0. It may be seen that zero of \Gcd(Q 2 ) \ predicted in ref. || is 
shifted to the lower values of Q 2 and the height of the secondary maximum is approximately 
the same as in our calculations. 

Here we restricted ourselves only to the discussion of the deuteron charge form factor 
Gcd- Even in this case we omitted such interesting questions as an analytical representation 
of relativistic corrections in different orders in (v/c) 2 , the new representation for realistic 
deuteron wave functions, the role of relativistic rotation of nucleon spins and orbital mo- 
mentum / = 2 in the deuteron, the problem of extraction, using the present approach, of 
Geu{Q 2 ) for ultralow values of Q 2 from experimental data on elastic ed-scattering, and con- 
tributions from meson-exchange currents. It would also be interesting to perform a detailed 
comparison of the present approach with other relativistic approaches to the description of 
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deuteron structure. 

All these questions, as well as the calculations of the deuteron magnetic and quadrupole 
form factors will be discussed in forthcoming publications. 
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APPENDIX A: KINEMATIC VARIABLES. 

By definition s is the invariant mass of n — p system squared: 

S = (Pn + Pp)l ■ 

In laboratory (LS) and center-of-mass (CMS) systems we have 

s = AM 2 + 2E = AM 2 + Ap 2 , 

where E is the nucleon energy in LS and p is the magnitude of the nucleon 3-momentum in 
CMS. 

Q 2 is the magnitude of the 4-momentum transfer squared: 

q 2 = = -t > o . 

APPENDIX B: JOST MATRICES B,B. 

The formulae for pairs (S, B) and (S, B) have the most convenient form in the p-plane: 

S(p)B+(p) = B_(p) , - oo < p < oo , 
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where S = S[5(p) , i](p) , e(p)}, see eq.©. Let us introduce two new matrices S and B: 
B±(p) =R(Tp)B±(p), 



t 



Rip) = I 



2ia 



(p+ia){l+p 2 ) 



\ 



(a 2 = Me). 



(Bl) 



Now the equation for B has the form 

S(p)B+(p) = B-(p), -o° <P<oo , 

S(p) = R(p)S(p)R-\-p) = S[6,i},e\. 
The last equation defines the reduced phase shifts 8, e, fj as functions of input experi 
mental phase shifts S,s,rj. 

The solution of eq.(BT) was found in ref. [rT| in the form of series 

oo 

B±(p) = B ±t0 (j>).[I+ y £B±, m (p)], 



m=l 



where 



( 



^ 1 ( p ) e T%) 

V2 (p)e T '' n{p) 



\ 



5(p')dp' 



1 r c 
<pi(p) = exp[ V.P. 

71 J-oo p — p ' 

f s r 1 T/P f°° WW\ 

(p 2 {p) = exp[ V.P. / 

7T J—\ 



-oo p — p 



1 



B±,m(p) = ~ 



dp' 



lixi J-oo p — p' ±iO 

In eg. (|B2|) for odd n 



J2G n (p')B +fi (p')[B + , m _ n (p>)} 



All-<5„ 



and for even n 



G n (p) = H-l)*? • ± ■ [2e(p))' 



G„(p)=i(-l)».i[2e(p)r 



\ 





e i(S+f/) q 

e 2fs 
e 2 ^ 



<5 mn is the Kroneker delta. 



(B2) 
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APPENDIX C: G^ L ' -MATRIX. 

In terms of invariant variables s, s', t and the nucleon EM form factors the matrix elements 
have the form: 

g(s, s', t)[gi(s, s', t)(cosai cos 0:2 — - sin ot\ sin a 2 ) • G S EN {Q 2 ) + 

o 

1 1 

+—g 2 {s,s',t) ■ (-sinaicosa 2 - cosaisina 2 ) • G S MN (Q 2 )] , 



9?{s, s',t) = 

g(s,s',t){gi(s, s',t)(-V2P 20 cosai sina 2 + -7=^21 sinai cosa 2 ) • G S EN 

v2 

— g 2 (s, s',t)(V2P 20 sina 1 cosa 2 H — 7=^21 cosai sina 2 )G s MN }, 



gl°(s, S ',t)=g° 2 ( S ',s,t), 



9? 'is, s',t) = 

1 2 

s',t){g 1 (s, s',t) [(-P21-P21 + 3^20^20) cos («i - "2) 
1 1 

+ (— -P22A2 + o-^20^2( 



+ 



0-P20-P20) cosai coso;2 + 
+ (y^(-P22-P 2 i - ^21^22) + 2^ 21 ^o ~ ^ 2 °- f 2i)) 



sin(o;i — a 2 ) — 



' G S EN - 2]^92(s, s',t) 



-7^(^22^20 + ^20^22) sin «i sin "2 

— ((-P2I-P22 ~ ^22-P 2 i) + -(-^20-^21 - -P2I-P20)) • 

cos(«i - a 2 ) - (^-^22-^22 + ^20^20) cosai sina 2 - 

-7^(^22^20 - P2qP 22 ) sinai cosa 2 + 

1 2 

+(-^21^21 + 3^20^20) sin(ai - a 2 )] • G^} , 
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where 

n ( Q q > f \ — gi(s,s',t)-(-t) _ i . i 

y{ ' ' ] v /(s-4M 2 )(s'-4M 2 ) [A(s, S ',t)]3/2 

#i(s, s',£) = s + s' -t, 

r n 1/2 

ff2 (s,s',t)= [(-l)(M 2 A( S , S ',t) + SS 't)] , 
A(s, s', i) = s 2 + s' 2 + t 2 - 2(ss' + st + s'i). 
p m are the Legendre polynomials, P/ m = Pi m {x) and P/ m = Pi m (x f ), where 

x(s s' t) = , V s ' (s '- S -V 

V ' ' ' A /(s'-4M 2 )A(s,s',i)' 

x'(s,s',t) — —x(s',s,t). 
The angles cki, «2 of the relativistic rotation of nucleon spins in deuteron are 

«i = arctan^ g f s ' s 3 _ , 

rvo = arrtan g2(s,5',f)( v ^+ v / 7+2M) 

2 d- 1 ^ 11 M( S + S '-()(^+V?+2M)+v^(4M2-() ■ 

r = Q 2 /AM 2 ; G S E ^ MN = \{G e ,m p + G E ,Mn) are the nucleon isoscalar charge and magnetic 
form factors. 
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FIG. 1. Separation of the unconnected part of the two-nucleon electromagnetic current. 
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FIG. 2. Neutron-proton phase shifts from ref. Q and their extrapolation. 
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FIG. 3. Relativistic deuteron charge form factor and its nonrelativistic limit. RDCFF with 
nonzero values of Geu = —^utGep is also shown. For comparison DCFF in relativistic approach 
of |8| and experimental results from MIT-Bates || are given, too. 
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FIG. 4. The various extrapolations of 3 Si-phase shifts to high energies. Curve 1 is the same as 
in Fig.2. The curves 2,3 are described in the text. 
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FIG. 5. RDCFF for the phase shifts from Fig.4. 
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